Fermi surface topology and de Hass-van Alphen orbits in Puln 3 and PuSn 3 compounds 
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Since the recent discovery of plutonium-based superconductors such as PuCoIns, systematic stud- 
ies of the electronic properties for plutonium compounds are providing insight into the itinerancy- 
localization crossover of Pu 5/ electrons. We are particularly interested in understanding the Fermi 
surface properties of Puln3 compound, which serves as the building block for the PuCoIns super- 
conductor. Motivated by the first observation of quantum oscillation and renewed interest in the 
de-Hass van- Alphen (dHvA) measurements on Puln3, we study the Fermi surface (FS) topology 
and the band dispersion in both paramagnetic and antiferromagnetic state of Puln3 , based on den- 
sity functional theory with generalized gradient approximation. We compare the results with its 
isostructural paramagnetic compound PuSn3. We show the detailed Fermi surfaces of compounds 
Puln 3 and PuSn 3 and conclude that the FS topology of an antiferromagnetic Puln 3 agrees better 
with dHvA measurements. In addition, we show the magnetization of the antiferromagnetic order 
can alter the angle dependence and values of the effective mass for the dHvA orbits. We point out 
the accurate determination of the magnetic order orientation with respect to the crystal orientation 
is crucial to advance the understanding of the electronic structures of the Puln3 compound. 

PACS numbers: 74.20.Mn, 71.27. +a, 74.70.Tx 



I. INTRODUCTION 



Actinide metals^l are strongly correlated electronic 
systems due to the narrow bandwidth of 5/ electrons. 
In addition, strong spin-orbit interaction within the 5/ 
electron systems comparable to other energy scales ren- 
ders the understanding of the electronic properties of 
the actinide metals more difficult. Itinerant-to-localized 
crossover of 5/ electrons that occurs near plutonium in 
the actinide series is one of the most challenging issues in 
condensed matter physics, partly because the dual char- 
acter (partially localized/delocalized) of these 5/ elec- 
trons is closely related to the abrupt atomic volume vari- 
ation between the a-Pu. and S-Pu metals. This change in 
bonding leads to a 25% larger volume in the S phase as 
opposed to a low-symmetry, monoclinic crystal structure 
a phase of Pu, along with a variety of unusual physical 
and mechanical properties. 3 To understand the 5/ elec- 
tron delocalization-localization crossover in elemental ac- 
tinide solids, it is very helpful to gain insight by studying 
its derivative compounds, which also show other emer- 
gent properties. One of the noticeable Pu-based com- 
pounds are 115 series. It has been discovered recently 
that superconductivity occurs in PuCoGas 4 (T c = 18.5 
K), PuRhGa 5 5 (T c = 8.7 K), and PuCoIn 5 6 (T c = 2.5 
K) series. The detailed pairing mechanism in Pu-115 
compounds are currently under intense investigation P 
We are particularly interested in the electronic properties 
of Puln3 and its isostructural partner PuSn3. The rea- 
sons are twofold. First, the PuCoIns structure consists 
of stacked C0L22 and Puln3 layers, therefore the under- 
standing of the electronic properties in Puln3 is relevant 
to uncovering the mystery of the superconductivity in 



PuCoIns. Second, the recent experimental capability to 
measure the Fermi surface topology of Puln 3 by dHvA 
effects can help narrow down the minimal effective theory 
for these complicated systems. 

In this paper, we present systematic studies of the 
electronic structures on Puln3 and PuSn3 within den- 
sity functional theory (DFT). We calculate the band 
dispersion, density of states (DOS), and the Fermi sur- 
face topology together with the identification of extreme 
orbits. For Puln3, there is recent experimental evi- 
dence 8 that the ground state of Puln3 is antiferromag- 
netic (AFM). However, the earlier dHvA measurement in 
this compound has been improperly interpreted in terms 
of DFT calculations based on a paramagnetic (PM) state 
of Puln3. 9 In the present work, we show the Fermi sur- 
face extreme orbits obtained for AFM Puln3 are in good 
agreement with the dHvA measurements, which places 
the theory and experiment comparison on a more con- 
sistent basis. Also with the known experimental fact of 
PuSn 3 being paramagnetic in the ground state, we pre- 
dict a quite different band structure and Fermi surface 
topology in PuSn3 as opposed to PM-PUL13. An ex- 
perimental verification of this prediction will provide a 
significant support for Landau's Fermi liquid theory and 
demonstrate the potential toward the detection of quan- 
tum oscillations in ^-phase of Pu. 

The organization of this paper is as follows. In Sec. II, 
we describe our theoretical method and briefly review the 
dHvA effects in itinerant electronic systems like metals. 
In Sec. Ill A, we discuss the band structures of the Puln3 
and PuSn3 compounds, respectively. In Sec. Ill B, DOS 
for the two compounds in concert to the band structures 
are discussed. In Sec. Ill C, we show the Fermi surface 
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FIG. 1: (Color online) (a) PM-PuM 3 (M=In,Sn) conventional 
unit cell. Red sphere: Pu, Grey sphere: M; (b) Top view 
of the anti-ferromagnetic Puln3 conventional unit cell; (c) 
Side view of the anti-ferromagnetic Puln3 conventional unit 
cell. The green arrows illustrates the commensurate AFM 
order. The tilt arrow just indicates the spontaneous symme- 
try breaking orientation for viewing convenience. The actual 
orientation for our studies can vary. 



topology and the corresponding dHvA orbits, which are 
relevant to available experimental studies. Discussions 
and conclusions are summarized in section IV. 



II. THEORETICAL METHOD 

Puln3 and PuSn3 compounds crystallize into cubic 
AuCu3-type structure at room temperature and have the 
actinide-actinide distance far above the Hill limit, 10 mak- 
ing the 5/-ligand hybridization the dominant mechanism 
for Pu 5/-electron derealization. The lattice constants 
are given by 4. 61 A for PUL23, 11 and 4.63A for PuSii3p21 
respectively. PuSn3 was reported as a paramagnet exper- 
imentally. 13 However, the Puln3 was originally assumed 9 
to be also paramagnetic but more recent experiments in- 
dicate it is in an AFM state below 14 K. 8 In this work, 
we consider both PM and AFM cases for the Puln3 com- 
pound. The conventional unit cell of PM-PUM3 (M=In, 
Sn) is shown in Fig. Il[a) with the ligand atoms M situ- 
ated at face centers of the cell. The magnetic unit cell 
for the AFM-Puln3 is enlarged because it is the com- 
mensurate AFM order obeying the discrete translational 
invariance in the crystal. Figure [lib) shows the top view 
of the magnetic unit cell accommodating the commen- 
surate AFM order. The side view is shown in Fig. Tic). 
Therefore, the lattices in this case can be considered as a 
tetragonal Bravais lattice decorated with the basis with 
Pu and In atoms. 

The first-principles calculations used here are based 
on the generalized gradient approximation with the PBE 
(96) functional!^ We use the scheme of the full-potential 
linear augmented plane wave basis plus local basis^ (FP- 
LAPW+lo), as implemented in the WIEN2K codeP 1 
The LAPW sphere radii R used for Pu, In, Sn atoms 
are 2.5 Bohr. The interstitial plane wavenumber cut-off 
kmaxR = 8-0 is chosen for the basis set. The semi-core 



states (Pu 6s6p, In 4p) are included with the valence elec- 
trons using local orbit als. The core states are treated at 
the fully relativistic level with spin-orbital interaction for 
all atoms. The spin-orbital interaction for semi-core and 
valence states are incorporated by a second variational 
procedure^ using the scalar relativistic eigenstates as 
basis, except that the so-called 6pi/2 relativistic local or- 
bitals 18 are used to account for the finite character of the 
wave function at the Pu nucleus. Dense Brillouin zone 
sampling is used for the calculation of the Fermi surface. 
We use a 31 x 31 x 31 k-mesh for the PM-Puln 3 and 
PuSn3 calculations. A 31 x 31 x 22 k-mesh is used for 
the AFM Puln3. Experimentally, a full determination of 
the magnetic structure for the AFM state of Puln3 is cur- 
rently not available. Therefore, we assume this stru cture 
follows that of its isostructural compound: CeIn3p^2Sl 
This choice is also supported by Fermi surface topology of 
PM-Puln3 as discussed later. The dHvA effect is observ- 
able in very clean metallic systems, typically in strong 
magnetic fields exceeding several tesla. On sweeping the 
magnetic field B, one observes oscillations in the magne- 
tization, which are periodic in an inverse magnetic field 
due to the fact that the number of occupied Landau lev- 
els changes with the magnetic field. The measurement 
of dHvA effect with varying magnetic field orientation 
is a powerful probe on Fermi surface topology in metal- 
lic and inter-metallic systems. The dHvA frequency F 
in MKS units is related to the extremal Fermi surface 
cross-sectional area A surrounded by extremal cyclotron 
orbits perpendicular to the magnetic field orientation!^ 



2ire 



(1) 



where e is the elementary charge. In addition, the effec- 
tive mass of electrons averaged around the cyclotron or- 
bits can be determined by the damping strength of dHvA 
measurements as a function of temperature. Our deter- 
mination of the dHvA orbits is based on the numerical 
algorithm implemented by Rourke and JulianP^J 



III. RESULTS 
A. Band structures 

Since PM-PUM3 compounds can be generated from the 
cubic unit cell in space, the corresponding BZ is cubic in 
reciprocal space. For simplicity, we label high-symmetry 
k points by the notation adopted by the tetragonal unit 
cell, as shown in Fig. [2j even though the BZ for the PM- 
PuM 3 is cubic. 

The band structure for both PM-PUL23 and PuSn3 
demonstrates two bands across the Fermi energy. The en- 
ergy dispersion along the high symmetry points is shown 
in Figs. [3] and [4j We observe the Fermi energy is crossed 
by two bands mainly of 5/ electron character, as shown 
later in the density of states (DOS) studies. The /-bands 
are far from flat due to the hybridization mainly from 




FIG. 2: (Color online) Tetragonal Brillouin zone (BZ). The 
symmetry points are labeled by capital letters. 





FIG. 3: (Color online) Band structure of the PM-Puln3 along 
the high symmetric points in BZ. 



FIG. 4: (Color online) Band structure of the PM-PuSn 3 along 
the high symmetric points in cubic BZ. Two conduction bands 
are across the Fermi level Ef- 



the 5p valence electron of ligand atoms. We calculate 
the band structure of PM-Puln 3 and PM-PuSn 3 by the 
density functional theory using the GGA approximation. 
We notice that the energy dispersion for PM-PuSn3 is 
very different, which indicates that the PM-Puln3 and 
PM-PuSn3 have very different Fermi surface topology, 
which is discussed in Sec. Ill C. 

For the AFM Puln3 with the chosen magnetic unit cell, 
the BZ in the reciprocal lattices becomes tetragonal. Let 
us first discuss the circumstances in which magnetic mo- 
ment orientation a for the AFM order is parallel to [001] 
direction defined by PM-Puln3 crystal unit cell. There 
are more atoms in the unit cell and therefore far more 
bands are accommodated in the corresponding tetrago- 
nal BZ. We observe four bands crossing the Fermi energy 
Ef as shown in Fig. [5] In this case, the bandwidth of 
5/ bands across Fermi surface is much narrower than 
the bandwidth in the PM case. The energy bands near 
the Fermi energy have dominant 5/ electron features, as 
shown later in the DOS, and have large effective band 
mass because of the flatness of the bands. 



B. Density of states 

By the smaller equilibrium volume of Puln3 in compar- 
ison with PuSn3, one may argue that PM-Puln3 should 
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FIG. 6: (Color online) DOS per atom and total DOS in PM- 
Puln3 compound. Ten times of the actual density of states 
for p orbitals at In site are plotted for better contrast. The 
black solid line shows the total DOS in the BZ. 



FIG. 5: (Color online) Band structure of the AFM-Puln 3 
along the high symmetric points in tetragonal BZ with the 
local spin polarization cr || [001] . 



have a more itinerant 5/ bands in comparison with PuSn3 
due to the smaller distance between nearest-neighboring 
Pu atoms and the In atoms. In addition, the atomic 5p 
levels from ligand atoms in Puln3 tend to be more de- 
localized due to smaller attraction from nucleon charges. 
However, those effects are not expected to be very re- 
markable in the DOS because of the closeness of the 
atomic properties for the In atom and Sn atom. 

To understand the hybridization between Pu 5/ and 
ligand valence states, we perform the DOS calculation 
for these compounds. In Fig. [6| the DOS for the PM- 
Puln3 is shown. We observe a nonzero total DOS near 
the Fermi level. The peak near the Fermi energy shows a 
hybridization of Pu 5/ with the DOS of 5p valence elec- 
trons from the In atom. This indicates that the metallic 
behavior is inherent from the hybridization of 5/ local- 
ized electrons at Pu with bp itinerant electrons at In. The 
two strong peaks away from the Fermi surface are due to 
the strong spin-orbital coupling for the localized / elec- 
trons, which is weakly hybridized with In 5p orbitals. The 
lower peak is mainly due to j = 5/2 sub-bands and the 
higher peak is mainly due to j = 7/2 sub-bands, which 
has been studied previously by relativistic Stoner the- 
ory)^ The 5/ electron occupancy at Pu atomic sphere is 
around 4.8. This suggests that the mid-peak near Fermi 



energy originates from the higher j = 5/2 sub-bands in 
the presence of crystal-field splitting. The three-peak 
features are similar also in PuSn3 except the j = 5/2 
peak are broadened by the hybridization from Sn p or- 
bitals. Fig. [7] presents the DOS for PuSn3. We notice 
the j = 5/2 sub-bands are still sharply peaked but with 
a slightly larger 5/-DOS at Pu site as shown by the solid- 
red line as suggested by a slightly smaller hybridization. 
As will be mentioned later, PM-PUL23 and PM-PuSn3 
have very different Fermi surface topology despite of the 
similarity in DOS. 

We have also calculated the DOS for the AFM-PUL23, 
as shown in Figs. 1(b) and 1(c), the lattice can be divided 
into two sub-lattices. The calculated spin projected par- 
tial DOS at different atomic sites is shown in Fig. [8] We 
observe that the Pul atom has larger 5/ occupancy in 
the up-spin channel (see corresponding partial DOS as 
denoted by "blue" solid line). However, the Pu2 atom 
takes larger 5/ occupancy in the down-spin channel (see 
corresponding partial DOS as denoted by "blue" dotted 
line). This indicates that Pul and Pu2 atoms have mag- 
netic moments due to 5/ electrons, which are equal in 
magnitude but antiparallel in orientation. We do not ob- 
serve net magnetic moment at In atoms according to the 
equal spin-projected DOS for the Inl atom indicated by 
the solid-dashed curves. The AFM order has the effect of 
generating a gap, above the Fermi energy, between two 
major peaks belonging to up and down spin Pu 5/-DOS. 
By our first-principle calculation, we found the orbital 
moment at Pul atom is around — 2.1/i# (jib '• Bohr mag- 
neton) but changes sign (-2.1/i#) for Pu2 atom. The spin 
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FIG. 7: (Color online) DOS per atom and total DOS in a 
PM-PuSn3 compound. Ten times of the actual DOS for p 
orbitals at Sn site are plotted for better contrast. The black 
solid line shows the total DOS in the BZ. 



moment at Pul and Pu2 atoms are given by 4.7/i# and 
— 4.7/xb respectively. We observe the presence of mag- 
netic order causes the mixing of j — 5/2 and j — 7/2 
and turn the local dip structure near the Fermi surface 
to a peak structure. In addition, two major peaks in 
DOS appear away from Fermi surface with large energy 
separation between maxima. 



C. Fermi surface topology and dHvA orbits 

In this section, we show the Fermi surface topology and 
calculate the dHvA frequencies as a function of magnetic 
field orientation for compounds PuSn3 and Puln3. By 
the band structure calculation, Fermi surface topology 
is determined by the isosurface calculated by the eigen- 
states of Kohn-Sham equation. 



PM-Pulri3 compound 

The Fermi surfaces for PM-Puln3 from two conduc- 
tion bands are shown in Figs. [9] and 10 The ordering 
of bands are based on the energies in increasing order 
in the band dispersion (Fig. |3|. There are two bands 
across the Fermi level. The characters of the dHvA or- 
bits are tabulated in Tables I and II respectively. We 
present the dHvA orbits of interest by "yellow" orbits 
or the label "orbit 1" for later discussion. This is high- 
lighted in the final two columns in Tables I and II. For 
Band-1, when the magnetic field is chosen along [111] 



FIG. 8: (Color online) Spin-projected DOS for the AFM- 
Puln3 compound. Pul represents the Pu atom with a positive 
local magnetic moment. Pu2 represents the Pu atom with a 
negative local magnetic moment. It is noted that Pul atom 
and Pu2 atom are the nearest neighbors in the AFM unit cell. 
The Inl atom indicates the indium atom located between the 
Pul and Pu2 atom in the unit cell (Figs. 1(b) and 1(c)). The 
color "blue" indicates the up-spin projected DOS along the 
chosen quantization -\-z axis in Bloch sphere. The color "red" 
indicates the down-spin projected DOS along the — z axis in 
Bloch sphere. Ten times of the actual spin-projected DOS of 
p orbitals at the Inl site are plotted for better contrast. 



direction, four dHvA orbits are identified in the BZ as 
listed in Table I. Three of the orbits are electron-like and 
one is hole-like. The dHvA orbit labeled by 1 originates 
from the Fermi surface near the center of BZ and has the 
dHvA frequency F = 0.9627 kiloTesla(kT). The location 
of the orbit 1 is shown in Fig. [9] explicitly. The other 
orbits with much larger dHvA frequencies are not our 
main focus since the dHvA frequencies are far from what 
have been observed experiment allyP For Band-2, only 
an electron-like (yellow) orbit with the dHvA frequency 
F = 3.9206 kT near the corner of the BZ is found, which 



is shown in Fig. 10 As far as the field angular dependence 
is concerned, we focus on the orbits which are connected 
with the same Fermi surface surrounded by the yellow 
orbits from Band-1 and Band-2 at the field orientation 
[111]. We observe the continuation of the orbits from the 
same Fermi surface spanned the whole angle range from 
[100] to [111] field orientation. In addition, the calcu- 
lated dHvA frequencies are not close to 2 kT as shown in 
Fig.[TT[a). Those facts contradict the observed dHvA fre- 
quency (around 2 kT) which only survives a limited span 
of field angles P This indicates that the assumption of 
the PM state in Puln3 is problematic for the interpreta- 
tion of the electronic properties of the Puln3 compound. 




FIG. 9: (Color online) Fermi surface from Band-1 in PM- 
Puln3 compound. 



From the Fermi surface topology for PM-Puln3 in Figs. 9 
and 10, we notice the nesting between the small Fermi 
surface at the T point from Band-1 and the four Fermi 
surfaces at the corners of BZ from Band-2. The nesting 
wave vector between the Fermi surfaces is approximately 
(7r/a, 7r/a,7r/a). Therefore, magnetic instability for com- 
mensurate AFM order should indeed be favored. In the 
following, we study commensurate AFM order in Puln3. 

TABLE I: dHvA orbits from Band-1: PM-Puln 3 , B || [111] 
T(kT) m*(m e ) Type Number Color Label 



16.149 
13.209 
10.818 
0.9627 



6.4320 
4.9218 

6.02 
2.4587 



TABLE II: The dHvA orbit from Band-2: PM-Puln 3 , B 

[111] 

F(kT) m*(m e ) Type Number of orbits Color Label 



3.9206 3.3774 e 



yellow 1 



AFM-Pulri3 compound 

From the band structure calculations as shown ear- 
lier in Fig. [5j there are 4 bands across the Fermi level 
for AFM Puln 3 . The labeling of bands still follow the 
band energies (Band-1, Band-2, etc.). The correspond- 




FIG. 10: (Color online) Fermi surface from Band-2 in PM- 
Puln 3 compound. 



ing Fermi surfaces and dHvA orbits (when the external 
magnetic field is orientated along [111] direction) are il- 
lustrated in Figs. [l2][T5l respectively. We are interested in 
dHvA orbits with dHvA frequencies around 2 kT, which 
are labeled as orbit 1 for Band-2 and Band-3 in Tables 
IV and V. The information of other dHvA orbits due 
to Band-1 and Band-4 are listed in Tables III and VI. 
We provide the Fermi surface topology for Band-1 and 



Band-4 in Figs. 12 and 15 for reference. 

For Band-2, the yellow orbits with the dHvA frequency 
of 1.3908 kT around the Fermi surface at the corner of the 
BZ are illustrated in Fig. 13 When the field orientation 



varies, one expect dHvA orbits at different angle are con- 
tinued orbits from the same Fermi surface. Therefore, we 
expect the dHvA orbits inherited from the same Fermi 
surface span the whole range of the angle as shown in 
Fig. 16 The dHvA frequencies are lowest around [111] 
orientation (angle= 54.7356°) and the value is around 
1.5 kT. The effective band mass is around 5m e where 
m e is the free electron mass. For band-3, in Fig. [l4j 
the dHvA frequencies around 2.3 kT spanned through 
a limited range of angles around [111] as can be under- 
stood by the shape of the corresponding Fermi surface. 
The dHvA frequencies have a local maxima around [111] 
orientation instead. The effective mass is above 10m e 
. Even though we have a good agreement with the ex- 
periment on the dHvA frequencies with a limited angle 
span, the predicted dHvA frequencies have a local max- 
ima instead of a minimum as observed experimentally. 
There is uncertainty in assuming magnetization orien- 
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FIG. 11: (Color online) Angular dependence of dHvA fre- 
quencies and the average effective mass of the yellow orbits in 
PM-Puln3 compound. is the angle spanned from the [100] 
to [011] orientation at the same plane. 




tation of AFM order along [001] orientation. We expect 
the orientation of the magnetic order will lead to qualita- 
tive effects on the Fermi surface topology because of the 
strong spin-orbital coupling effects in Puln3 compounds. 
Therefore, we have also considered an easy axis of mag- 
netization along [111] direction. In this case, we predict 
three 5/ bands across the Fermi level in this case. We 
show the Fermi surface topology corresponding to three 
bands across Fermi level in Figs. [l7}fT9l The details of 
the dHvA orbits along [111] orientation are shown in Ta- 
bles VII, VIII, and IX. We identify an orbit with dHvA 
frequency F = 2.125 kT marked by a yellow trace in the 
Fermi surface around T point from Band-2 in Fig.[l8] By 
the calculated angular dependence of the dHvA orbits for 
Band-2 in Fig. 20 'a), we identify the dHvA orbits are lim- 
ited to certain angle range around [111] orientation. In 
addition, the effective mass has a local minima around 
[111] orientation as suggested by experimentP However, 
the effective mass is only about 2m e . For the dual nature 
of the 5/ electrons, the temporal quantum fluctuations 
can renormalize the effective band mass dramatically. We 
do not expect density-functional theory as presented in 
this paper can capture this effect. However, the Fermi 
surface topology and thus the predicted dHvA frequen- 
cies should still be robust by Landau's Fermi liquid the- 



FIG. 12: (Color online) Fermi surface from Band-1 in AFM- 
Puln3 compounds with the local spin polarization cr || [001]. 




PM-PuSri3 compound 

So far, dHvA effects have not yet been successfully ob- 
served in PuSn3 compounds. By our calculation, we find 
multiple dHvA orbits in a narrow range of frequencies at 



FIG. 13: (Color online) Fermi surface from Band-2 in AFM- 
Puln3 compounds with the local spin polarization cr || [001]. 



TABLE III: The dHvA orbit from Band-1: AFM-Puln 3 , B 
[111], a || [001] 



F (kT) m*(m e ) Type Number of orbits 




0.4424 
0.4412 



2.1654 
2.3620 



TABLE IV: dHvA orbits from Band-2: 
[111], a || [001] 



AFM-Pulns, B 



F (kT) m*(m e ) Type Number of orbits Color Label 



1.3908 6.2948 

1.006 3.7899 

0.5189 2.6069 

0.3335 2.2621 



yellow 



TABLE V: dHvA orbits from Band-3: 
[111], a || [001] 



AFM-Pulns, B 



F (kT) ra*(ra e ) Type Number Color Label 



FIG. 14: (Color online) Fermi surface from Band-3 in AFM- 
Pulns compounds with the local spin polarization cr || [001]. 



2.3615 
0.7796 
0.6833 
0.577 
0.2279 
0.0386 



13.4852 
1.8097 
3.0840 
11.017 
2.2562 
0.8525 



yellow 1 



TABLE VI: The dHvA orbit from Band-4: AFM-Puln 3 , B 

[111], a || [001] 

F (kT) m*{m e ) Type Number 



0.5277 1.0592 
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FIG. 15: (Color online) Fermi surface from Band-4 in AFM- 
Pulns compounds with the local spin polarization cr || [001]. 



FIG. 16: (Color online) Angular dependence of dHvA frequen- 
cies and effective mass of chosen yellow orbits for AFM-Puln3 
compound. 



^ * ^ 




FIG. 17: (Color online) Fermi surface from Band-1 in AFM- 
Puln3 compounds with the spin polarization <x || [111]. 



FIG. 19: (Color online) Fermi surface from Band-3 in AFM- 
Puln3 compounds with the spin polarization <x || [111]. 




TABLE VII: The dHvA 
B || a || [111] 


orbit 


from Band-1: AFM-Puln 3 


F (kT) ra*(ra e ) Type 


Number of orbits 


Color Label 


3.2569 5.8858 
1.7165 1.9391 
0.9876 3.4407 
0.9697 2.8944 
0.8421 3.5473 


h 
h 
h 
h 
h 




1 
1 
2 
2 
2 


yellow 1 


TABLE VIII: dHvA orbits from Band-2: , 

a || [111] 


\FM-Puln 3 , B | 


F (kT) ra*(ra e ) Type 


Number of orbits 


Color Label 


3.7022 6.2443 
2.1250 1.9877 
1.2687 2.8742 
1.1595 1.5520 
1.0834 1.7674 


e 
e 
e 
e 
e 




1 
1 
2 
1 
2 


yellow 1 



TABLE IX: dHvA orbits from Band-3: AFM-Puln 3 , B 

[111] 

F (kT) ra*(ra e ) Type Number 



FIG. 18: (Color online) Fermi surface from Band-2 in AFM- 
Puln 3 compounds with the spin polarization <x || [111]. 



0.3359 2.3248 

0.3358 2.4819 

0.3056 2.1990 

0.1935 1.7949 



10 



(a) 



Local spin polarization || [111] 



2.5 



(>- 



1.5 L 
20 



- Band 1 : orbit 1 
Band 2 : orbit 1 



-e- 



(b) 



-e- 

30 40 50 

Angle(degree) 




60 



70 




40 50 

Angle(degree) 

FIG. 20: (Color online) Angular dependence of dHvA fre- 
quencies and effective mass of the chosen yellow orbits for 
AFM-Puln3 compound. 



the field orientation [111]. This may cause complications 
on observing clean signals. However, we do find a sim- 
ple Fermi surface topology at the field orientation [001], 
which can be easier to observe and thus is the case we 
focus on in this section. With the PuSn3 band struc- 
ture in Fig. [4j we observe two energy bands crossing the 
Fermi energy Ep, which only have similar dispersion as 
PM-Puln 3 near the point Y = (0, 0, 0) in BZ. The cutting 
points at the Fermi energy in Fig. [4] are the projection of 
the three dimensional Fermi surface plot along the line 
between symmetry points in BZ. We show the Fermi sur- 
faces for those two bands in Figs. [21] and [22] respectively. 
Our calculation shows one hole-like (yellow) orbit (with 
F = 3.273 kT ) near the Y point in the BZ, which agrees 
with the band dispersion at Band-1 near the Y point. 
For Band-2, the four dHvA orbits are identified and all 
of them are hole-like at the field orientation [001] (See 
Table VIII) , denoted by colored closed curves in Fig. [22] 
The maximum dHvA frequency is given by the white cy- 
clotron orbit with F = 1.312 kT near the M point in BZ. 
By observing the Fermi points and the band dispersion 
nearby for Band-2 in Fig. [4j one can see the agreement 
with the Fermi surface. In addition, from the Fermi sur- 
face topology, one can see there is no clear evidence for a 
dominant Fermi surface nesting. This might explain the 
absence of magnetism in PuSnsP^ 

In Fig. [23] we show the field angle dependence of dHvA 
frequencies for the yellow dHvA orbit in Band-1 and the 
white dHvA orbit in Band-2. The magnetic field varies 
from the crystal orientation [100] to [001] in the PuSn3 
compound. We observe that the dHvA frequencies are 
symmetrical with respect to the orientation at the angle 
45°. The dHvA frequency for Band-1 is insensitive to 
the orientation because of the corresponding spherical 




FIG. 21: (Color online) Fermi surface from Band-1 in PM- 
PuSn3 compound. 



TABLE X: The dHvA orbit from Band-1: PM-PuSn 3 , B 
[001] 



F (kT) m*(m e ) Type Number Color 
3.273 L5976 h 1 yellow 



TABLE XI: dHvA orbits from Band-2: PM-PuSn 3 , B || [001] 



F (kT) 


m*(m e ) Type 


Number of orbits Color 


0.5550 


0.4749 h 


2 yellow 


0.6873 


0.7955 h 


4 blue 


1.230 


1.6842 h 


2 green 


1.312 


0.8721 h 


1 white 



Fermi surface near the Y point and is well separated in 
frequencies from the orbits due to Band-2. However, the 
Fermi surface circled by the white orbit in Band-2 is not 
spherical (shown in Fig. 22); therefore, a strong field 
angle dependence due to the Fermi surface is obtained. 



IV. CONCLUSIONS AND DISCUSSIONS 

We have studied the electronic properties for com- 
pounds PuSn3 and Puln3 within density functional the- 
ory under the GGA approximation. Our calculations for 
AFM-Puln3, with magnetization along [001] direction, 
show that the Fermi surfaces with heavy effective band 
mass produce the dHvA frequencies near 2 kT, in quan- 
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titative agreement with experiments but the angular de- 
pendence of dHvA frequencies around [111] orientation 
cannot be reconciled with experiment. However, when 
the magnetization is parallel along [111] direction, we 
can identify new orbits showing good agreement with ex- 
perimental measurement in dHvA frequencies. We point 
out the predicted effective band mass values can be fur- 
ther renormalized due to temporal quantum fluctuations 
which cannot be captured properly with the density- 
functional theory. We suggest experimental investiga- 
tions of the effective band mass of the dHvA orbits and 
the magnetic order are crucial to resolves those issues. 
Finally, we have also predicted the band structure and 
dHvA frequencies for PuSn3, which is known to be para- 
magnetic, for guiding the search of dHvA effects in the 
PuSn3 inter-metallic compounds. Experimental observa- 
tion of these frequencies in the compound will help the 
search for quantum oscillations in £-phase of Pu. We 
do not expect the predicted values for the effective mass 
is reliable with density-functional theory ignoring strong 
quantum temporal fluctuations but the predicted Fermi 
surface should be robust with respect to the given mag- 
netic order. 



FIG. 22: (Color online) Fermi surface from Band-2 in PM- 
PuSn3 compound. 
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FIG. 23: (Color online) Angular dependence of dHvA frequen- 
cies in PM-PuSn3 compound. The dHvA orbit of Band-1 is 
inherent from the same Fermi surfaces circled by the yellow 
orbit in Fig. 17. The dHvA orbit of Band-2 at different angle 
is inherent from the same Fermi surfaces circled by the green 
orbits in Fig. 18. 
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